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AND COOLING OF A BAFFLED CYLINDER BARREL 
By Arthur W. Goldstein and Herman H. Ellerbrock, Jr. 



SUMMARY 

Theoretical investigations have shown that, because air is com- 
pressible, the pressure-drop requirements for cooling an air-cooled 
engine will be much greater at high altitudes and high speeds than 
at sea level and low speeds. Tests were conducted by the NACA to 
obtain some experimental confirmation of the effect of air compressi- 
bility on cooling and pressure loss of a baffled cylinder barrel and 
to evaluate various methods of analysis. The results reported in the 
present paper are regarded as preliminary to tests on single -cylinder 
and multicylinder engines. Tests were conducted over a wide range of 
air flows and density altitudes. 

The results indicate that, for a given air weight flow, the 
reduced pressure drop based on average cooling-air density App av /p s , 
which has been used to correlate heat-transfer cooling data, is not 
constant for different air densities (Ap, cooling-air pressure drop 
across the engines Pav/Ps* .ratio of average cooling-air density to 
density at sea level). Engine-cooling variables should therefore not 
be plotted against pressure drop. From the present tests a correlating 
variable for heat-transfer data is shown to be the air weight flow; 
the reduced pressure drop is not suitable for this purpose. An analysis 
based on the assumption of uniform flow is shown to be satisfactory 
for estimating the effect of compressibility on data obtained in these 
tests. A simpler empirical method in which compressibility and heating 
effects, can be estimated was found for correlating the test data on 
pressure loss. 



INTRODUCTION 

Some investigators have heretofore correlated cooling data for 
air-cooled engines with the cooling-air pressure drop (references 1 
to k) . The effect of compressibility was taken into account by using 
the product of the pressure drop and the average of the air densities 



\ 



2 MCA AERNo. E4G20 

at the front and the rear of the cylinder as the correlating vari- 
able rather than the pressure drop alone. With the high rate of heat 
exchange and the high air velocity "between the fins required, for 
effective cooling at high altitudes, however, a largo air-density 
change will result. This change in air density is attended by an 
increase in velocity, and an additional pressure loss will he incurred 
at the baffle exit where this momentum will he lost. Estimates of 
the increase in pressure loss caused "by air compressibility in engine- 
cooling systems with baffled cylinders wore made in references 5 to 9. 
The analyses of references 5 and 9 were based on the assumptions of 
one -dimensional, gas dynamics, but no experimental data were available 
to support these analyses.' 

An investigation was begun by the NACA to obtain experimental 
confirmation of the effect of compressibility and rate of heat transfer 
on pressure less and cooling of a baffled section of a cylinder barrel 
and to evaluate the various methods of correlating these data. The 
tests covered a range of simulated density altitudes from 4000 to 
33,000 feet, of velocities between the fins from a Mach number of 0.05 
to near sonic values, and of heat inputs from 0 to 500 3tu per hour 
per square inch cf cylinder- wall surface. This investigation was con- 
ducted at Langley Memorial Aeror=auti.cal. laboratory— Langley-Fxeldy Va,, 
during 1942. 

Acknowledgment is made to Mr. Frank E. Marble of the Supercharger 
and Airflow Research Division, Aircraft Engine Research Laboratory of 
the NACA for his suggestion that the Prandtl-Glauert compressibility 
factor be used to correlate the results of the present tests. 



ANALYSIS 

One-Dimensional Gas Dynamics 

The analysis of the flow around a baffled cylinder is based on 
the assumptions of one- dimensional gas dynamics. The equations for 
the pressure and density changes through a baffled- cylinder system 
(fig. l) are developed by the following analysis, which is similar to 
that of reference 5. 

1. The heat transfer between stations 1 and 2 and between sta- 
tions 2 and 4 is estimated. 

2. No loss in total pressure is assumed from station 1 to sta- 
tion 2. 

3. The pressure and the density at station 2 are then calculated 
from the heat-transfer estimate, the stagnation pressure, and the 
mass flow. 
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U. The analysis of the flow between stations 2 and 3 is based 
on the assumption of uniform flow of a compressible fluid with fric- 
tion in a straight duct. 

The loss at the baffle exit is computed from the momentum 
equations for uniform flow at. stations 3 a ™i U. 

Change in gas state between station 1 and station 2. - If the 
• heat picked up by the air between station 1 and station 2 and the 
weight of air flowing through the baffle are known, 

% = wc p (r 2t - T lt ) (1) 

where 

H n rate of heat transfer to air between stations 1 and 2, 3tu per 
second (The method of obtaining is given in appendix A.) 

W weight of air flowing through baffle, pounds per second 

Cp specific heat of -air at constant pressure, Btu per pound per °F 

To stagnation temperature at station 2, °F absolute 

stagnation temperature at station 1, °F absolute 

k complete list of the symbols used is presented in appendix B. 

Stagnation or total temperature and pressure as used in this 
report indicate gas properties that would obtain if the kinetic energy 
of the moving gas were isentropically converted into enthalpy. The 
temperature T^ may be considered equal to the static temperature 

at station 1 because the velocity at station 1 is negligible. If Tn 
1, and ft, are known, the total temperature at station 2 may be 
obtained from equation (l). 

The assumption is made that there is no loss in total pressure 
from station 1 to station 2 and that the stagnation state of the air 
at station 2 is kno^n. Then, from the relations 

W 

and 

.do 

c 



P2V 
YP 2 



(2) 
(3) 



k 
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and by means of the relations for isentropic change, the true-stream 
density and pressure may be eliminated and the stagnation density 
and pressure inserted to give 



(p V 2 ) 2 YMp 

- ■ (U) 



■ 



1 + Lr-i M 2 2 



where 

2 

g acceleration of gravity, feet per second 

A 2 cross-sectional area of free-flow soace at station 2, square 

feet 

p 2 density of air at station 2, slugs per cubic foot (based on p? 
and T 2 ) 

V 2 velocity of air in fin passage at station 2, feet per second 

M 2 Mach number at station 2 

y ratio of specific heats for air (L.y)k7) 

P 2 static pressure of air at station 2, pounds per square foot 

Pp stagnation air pressure at station 2, pounds per square foot 

pg._ stagnation air density at station 2 (computed from p 2 and 
" t slugs per cubic foot " * 

T 2 true air-stream temperature of cooling air at station 2, 
°F absolute 

When (9^2^ 2 / P 2 t p 2 t is known, M 2 can be determined from equa- 
tion (h)i the ratios P2/P2 t ancl p 2/ p 2 t can te determined from the 

isentropic-change relations. A chart based on equation (I4) and the 
isentropic relations is given in figure 2. 

Change in gas state between sta ti on 2 and station 5> - The pres- 
sure drop through the baffle from station 2 "to station 3 is given by 
the following equation, which has been modified from the corresponding 
equation in reference £ by inserting the value of C Cf in terms of 

G 0f :L ' 
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P2 - ?3 _ C ° f i ^ p 2 N x 0 f" p 2 



q 2 

where 



q 2 



dynamic pressure at station 2, ( -Po^o J * pounds ner square 
foot * \2 * * / 

friction-drag coefficient between stations 2 and 3 based on 
friction-drag coefficient between station 2 and station 3 



based on average of and q*, Tc^, ~ (Df/kpJ f- — 

|^ *i 1 + ^3/ 

Djf friction drag from station 2 to station 3i pounds 

q? dynamic pressure at station 3* (ip?V^) , pounds per square 
9 foot \ 2 J ? / 

p^ static pressure at station 3* pounds per square foot 

density of air at station 3> slugs per cubic foot 

velocity of air at station 3* feet per second 

If the value for the local friction-drag coefficient C^, is 

assumed constant for all elements of the path, Df would actually be 
obtained by a process of integration. The value of Cr^. should 

therefore be calculated from Df by means of some integrated mean 
IP 

value for ^ va ^ ue calculated from the arithmetical average 

of and q^, however, is used as an approximation. If this 

approximation is good, Cr> P thus calculated should be independent 

"J-i 

of inlet-density variations for fixed values of Reynolds number and 
equal to the value that would be obtained with an incompressible 
fluid. By means of the continuity equation, the relationship between 
Cjjj and Cq^. can be established as 

1 ■ ^-^-d. <6) 
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The density/ ratio required for the solution of equation (5) can 

be calculated from the following equation, which is modified, from 

the corresponding aquation in reference $ by inserting the value for 

C;> in terms of Ctv>. to give 
x x 

- \ 

2 Or 



p 3 



=1 



n 2 Y-l „ 2 
+ __± Mo 



1 

1 "V.P3, 



*3 



- 1 



(7) 



-'mere 



H 2 rate of heat transfer to air between stations 2 and 5, Btu per 

second (The method of determining H 2 is given in appendix A.) 

A convenient method of determining P2/P3 uses the variables 



Ho 



T« = 



o ;,i 2 



(8) 



P 1 = p 2 /pj 



(9) 



M 2 . 



p' - 1 



Tpi 



L « 



-i (1 + pi) + (p< - 1) 



(10) 



p' - 1 



(11) 



When use is made of these definitions, equation (?) becomes 



M..2 = M 2 + T , „ 



-"a 



D 



(12) 



Because M a * and M> 0 ^ are f mictions only of p f and Cj)f. f for 
each value of C 2f . the variable M a 2 may be plotted against Iv^ 2 
ana p ! may^be plotted against M a by use of the definitions for 



2 



and 



m, 2 

,4 b 



These plots are shown in figure J. The value M a 2 
must, however, be the ordinate satisfying the linear equation (12) 
with 



M 2 - as the y intercept, as the abscissa, and T f as 

the slope. The intersection of this line with the curve of M* 2 
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9 9 

against for the given C Df determines M a - . From the curve 

of M a 2 against p 1 for the given C Df , the value of p ! can be 
obtained • 

C hange in gas state between station 3 and station U « - The 
momentum equation characterizing the pressure loss between station 3 
and station h is derived for the baffle with a tailpiece. This loss 
as given in references £ a * id 9 was for a baffle without a tailpiece, 
where station k was in a section of very large area with very low 
velocity. For the present case the momentum equation is 

P 3 A 3 sin 63 - vkHi. + fi d3 x ~ D 3 " I (V U ~ V 3 sin e 3 } (l3) 

where 

6^ angle between radii of cylinder to cylinder rear and to station 3, 
degrees (See fig. 1.) 

Pi static pressure at station h, pounds per square foot 

A3 cross-sectional area at station 3, square feet 

A^- cross-sectional area of baffle exit at station I4, square feet 

p static pressure, pounds per square foot 

dS x projection of any element of cylinder-wall surface and of curved 
part of baffle surface back of station J on plane of A^, 
square feet 

D3 component of drag force normal tQ plane of and effective 
between stations 3 a nd h$ pounds 

Vjj velocity of air at station k, feet per second 

The integral of pdS is taken over the entire surface bounding 
the fluid between stations 3 and U except at flow cross sections of 
stations 3 and It. 

It can be shown that 

JdS x = ki - A3 sin 83 (lU) 

The pressure gradients will be proportional to the kinetic energy at 
station 3. Therefore, if 
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a ? J <P - V (15, 



if the coefficient for friction drag in the baffle exit is 

C 5 = ■ \ (16) 
2 P3 v 3 A U 

and, if these symbols and the continuity equation are employed, the 
momentum equation may be expressed as 

!LL2i . c, - a, - 2 % sin fl, ♦ (17) 

In this equation the term 2(A^/A[ i ) sin 9^ indicates that the 

momentum loss paused by the fact that the stream at the rear of the 
cylinder is not directed straight back cannot be neglected. This 
direction of the stream, however, tends to increase the value of 
and the pressure recovery. - TITe Recovery cbsfi feient ^» therefore, 

in some measure compensates for the fact that G^ ^ 90°. In the case 
of Aj<<A[ i , all the terms become small and consequently ' p^ is 
approximately equal to pj^. 

For the present report the losses across the cylinder will be 
represented by equations (£) and (l?). In the actual practice of pre- 
dicting the pressure losses, once the coefficients of equation (17) 
are evaluated, the density change may be evaluated from equation (17) 
and the energy equation. 



Application of Prandtl-Glauert Compressibility Factor 

to Flow across a Baffled Cylinder 

The foregoing theory for flow of a compressible fluid is appli- 
cable only to the case of uniform velocity or one-dimensional flow. 

For two-dimensional flow, the Prandtl-Glauert factor J 1 - M G 2 is 
used to compute the effect of compressibility. The Prandtl-Glauert 
factor is strictly applicable in flow conditions quite different from 
those existing around, the baffled cylinder now being considered; there- 
fore, this name should not be used for this application of the factor 
but will bo used for convenience. The proper application of this 
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factor is to static-pressure variations in the flow field of a body 
causing small perturbation velocities in an infinite uniform flow 
field with a frictionless, compressible fluid and with ho heat 
transfer. In that case the difference in static pressure at any 
point in the flow field from what the static pressure would be with 
an incom pressible - fluid can be computed by the Prandtl-Glauert 

factor </l - M Q 2 in the equation (reference 10) 

_ P - Pq _ Pj - P 0 i 

c pm -IT " ~ ' (1<3) 



or 



where 



P o pressure of fluid that is characteristic of flow, pounds per 
square foot 

q 0 dynamic pressure of fluid that is characteristic of flow 

^2 P ° P ounds per. square foot 

M 0 Mach number that is characteristic of flew 

p i static pressure at same point as p with .incompressible fluid, 
pounds per square foot 



density of fluid that is characteristic of flow, slugs per 
cubic foot 



V 0 velocity of fluid that is characteristic of flow, feet per 
second 

Cp^ pressure-loss coefficient with compressible flow 

C n . pressure-loss coefficient with incompressible flow 
^i 

The factor - M 0 2 was used to reduce the .total-pressure-loss 
data to values that might be expected without compressibility effects. 
The application of the factor to the present data is to be regarded 
solely as an empirical method of correlation. No rational basis for 
the use of this factor is presented herein.. 
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The pressure drop used in place of p - p in equation (18) 
for the baffled cylinder is ^ - p^. The method of calculating 

q 0 and M 0 was empirically determined by finding the pressure and 
the temperature that would give the best correlation of the data 
with different densities and various rates of heat transfer. The 
factor Cp was then computed 

p it ~ P ^t r— — 5- , x 

C p . * -L- Jl^uJ (20) 

The total-pressure loss to be expected with an incompressible fluid 
Ap^ is then 

"For the same mass flow (pV) but with standard density p. C 7 , 

Pi 

will remain constant because it depends on only the Reynolds number. 
Then 



where 

Ap. loss in total pressure under standard density conditions, 
" s pounds per square foot 

p s standard density (at 29.92 in. Hg and 60° F), slugs per cubic 
foot 

From equations (20), (21), and (22) 



& p . A p t ^/l - M 0 2 (23) 
s p s 

where 

A.P t E P lt - V ht 



The Effect of Compressibility on Heat Transfer 

The effect of compressibility on. the heat-transfer coefficients 
is estimated by the effect calculated for a flat plate by an equation 
given in reference 11. If the temperature of the plate is assumed to 
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be 300 F (an average fin and cylinder temperature to be expected 
with a k<0° F rear spark-plug temperature), the Reynolds number is 
assumed to be J>0Q0, and the free-strearn air temperature is assumed 
to be -67° F (the temperature at high altitudes where compressibility 
effects in air-cooled engines may become critical), the effect of 
the free-stream Mach number M on the local heat-transfer coefficient 
h >: is given by 

h x S (1 + 0.0^6 M 2 ) 

where h is the local heat-transfer coefficient that would exist 
x i 

with an incompressible fluid. The effect will be slightly less than 
shown in the foregoing equation for two reasons: (l) the value h v 

i 

is proportional to the local skin-friction coefficient, which in ref- 
erence 12 is shown to decrease slightly with increase of Mach number; 

(2) the highest possible value of M 2 = 1 will not exist over the 
entire cylinder. In practical cases M 2 will not be unity at any 
point around the cylinder. Consequently, the effect of compressibility 
on heat-transfer coefficients for the cylinder can be expected, to be 
negligible . 



APPARATUS AND TESTS 
Test Setup 

The copper-plated steel barrel section used in the tests was 
ljjj inches long with a SjLinch bore. The fins on the cylinder were 

O.^O inch wide, O.O36 inch thick, and spaced 0.105 inch. Inside the 
cylinder was a grate to aid the pickup of heat and to reduce the tem- 
perature variation around the inside of the cylinder. A metal baffle 
with a 6-inch tailpiece was fitted around the cylinder (fig. l) . The 
tailpiece allowed the cross flow at the back of the cylinder to 
diminish sufficiently to permit more reliable pressure reading* . The 
cross-sectional area of the exit of the baffle was 1.6 times the free- 
flow area between the fins. The unit was placed in an asbestos-lined 
metal box (fig. h) and sealed at all edges with furnace cement to 
prevent air leakage. 

The source of heat was an oil burner with a capacity of 1 or 
2 gallons of oil per hour, depending on the burner nozzle used. A 
firebrick furnace provided the space for the combustion, to be com- 
pleted before the hot gases came into contact with the cylinder grate. 
An auxiliary blower was needed with a large nozzle to supply the 
necessary air for combustion. 
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The flow of cooling air was created by two compressors used as 
vacuum pumps and operated in series. Each pump was driven by an 
85-horsepower engine. In order to prevent surging and to smooth out 
the flow of cooling air, large tanks were placed upstream and down- 
stream of the test cylinder. Bleed valves were placed in front of 
and between the pumps to provide fine control of the air flow. A 
tank with thin-plate orifices in each end was placed upstream of the 
test cylinder to measure the quantity of cooling air. A throttle 
placed between the orifice tank and the upstream surge tank v/as used 
to control the pressure of the cooling air in front of the test cyl- 
inder. A diagram of the aj>paratus is shown in figure 



Instrumentation 

Surface temperatures of the cylinder were obtained at 2k points 
(fig. 6) by means of insulated 28-gage chromel-alumel thermocouples 
spot-welded to the steel. The cold junctions of the thermocouples 
were inserted in a sealed wooden bo:-:. The temperature in the box was 
obtained with an alcohol-in-glass thermometer j the thermocouple poten- 
tial was measuredjwijth_a i^t*:mta^meier. Xhermo_couples wBre_also_use_d 
to measure the orifice temperatures and the cooling-air inlet and 
outlet temperatures. The cooling -air outlet temperature was measured 
downstream of the baffle tailpiece in an expanded section; it was 
therefore unnecessary to correct the readings for air velocity. The 
outlet duct was lagged to prevent heat loss. The accuracy of the 
temperature measurements was within' ±1? F. 

Total and static pressures were obtained in the tailpiece and 
at two stations on the cylinder; a total pressure was obtained imme- 
diately in front of the cylinder (fig. 6). The total -pressure tubes 
of 0.030-inch-diameter steel tubing with a 0.006-inch-diameter hole 
in the side of each tube were inserted vertically around the cylinder. 
The static pressures around the cylinder were measured by means of 
vertical tubes inserted through holes tapped into the fins (fig. 7). 
Station 2 was at the baffle constriction and station 3 1 was somewhat 
behind the baffle-expansion point. Because of the location of the 
tube, the static-pressure readings for station 3 f were not used in 
the calculations. Also, in view of the small distance between 3 ! and 
3, the reading p^, was used for p^ • Conventional -type pressure 

tubes were not used between the fins because it was thought they 
would block too much of the channel. A rake of conventional total- 
pressure tubes was used in surveying the total pressure in the tail- 
piece, and a wall tap was used to measure the static pressure (fig. l). 
The pressures were read on vortical -water or mercury manometers, 
depending on the range of pressures being measured. 
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The very small pressure drops across the thin-plate orifices 
were measured with a micromanometer . 



Tests 

Tests were conducted over a range of air flows at several NACA 
standard density altitudes from U000 to 33,000 feet simulated at the 
front of the cylinder. At each altitude, data were obtained with 
and without heat transfer for Mach numbers ranging from low values 
to the highest values obtainable with the apparatus. The highest 
Mach numbers occurred at the baffle exit, but no accurate data were 
obtained at that station. In the following table are given the 
ranges of Mach numbers at the baffle entrance M2 and at the baffle 
exit M 0 ; the M 0 values were computed on the assumption of no 
friction up to that point and are therefore lower than the maximum 
Mach number of the flow system. 



Cylinder 


Density 










condition 


altitude 


M 2 




M 0 






(ft) 










Without 


14,000 


0.0^ to 0.68 


0.0£ to 


0.68 


heat 


15,000 


' .12 to 


.70 


.12 to 


.70 


transfer 


23,000 


.17 to 


.68 


.17 to 


.68 




33,000 


.18 to 


.60 


.18 to 


.60 


With 


u,ooo 


.16 to 


.73 


.18 to 


.77 


heat 


III, 000 


.16 to 


.67 


.17 to 


.69 


transfer 


2i|,000 


.21 to 


.69 


.22 to 


-7U 




32,000 


,29 to 


.56 


.31 to 


.59 



Pressures and temperatures were recorded only after the test 
cylinder had practically reached a state of thermal equilibrium. 
The maximum allowable rate of change of temperature was about 3° F 
in £ minutes for the maximum temperature before readings were taken. 
The cooling-air outlet temperature was read before and after reading 
the cylinder temperatures and an average of the two air temperatures 
was used. The two readings in no case differed by more than 2 per- 
cent of the temperature rise. Temperatures were read to within 0.[| O F 
and pressures to within 0.01 inch of mercury or water. The accuracy 
of the pressure readings was limited by the unsteadiness of the 
engines, which caused pressure fluctuations of as much as 0.1 inch of 
mercury . 

Summaries of the reduced test data and the derived quantities 
without and ytith heat transfer are given in tables 1 and 2, respec- 
tively. Cylinder temperature-distribution data are available upon 
request from the NACA. 
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COMPUTATIONS 

From the tests the following quantities were determined: 
W, H x + H 2 , p lt , Ti t , p 2 , p 2t , p 5t , p^, p^, T^, T^, T b , a, 83, 
A 2 , A3, A^, A b , arxd S T 
where 

A b area of cylinder wall at base of fins, square feet 

T b average wall temperature, °F absolute (arithmetic average of 
measured temperatures of cylinder at base of fins) 

T f average temperature of cooling surface (fin and barrel surfaces), 
av °F absolute 

S T total heat-transfer surface of cylinder, square feet 

From these data, the fictitious exit- density p and the coeffi- 
cients h, U, Cq^ t j'Cp., and -"€3 "were "determined as follows: 

Computation of h. - The average surface heat-transfer coeffi- 
cient h was calculated from the equation . 

h * S T ( T f av~ T lt) ' 
The average temperature Tf of the fin and barrel sui'face was obtained 

from the measured surface temperatures by averaging the temperatures 
with weighting factors proportional to the area elements in which each 
thermocouple was located. 

Computation of U. - The average wall heat-transfer coeffi- 
cient U was calculated from the following equation 

II-. + Hp> 

U = — (2W 

A b (Tb ~ Tl t ) 

Co mputation of fictitious exit density p . - The stagnation 
pressure at station 3 T , the mass flow, the exit stagnation tempera- 
ture, and the cross section at station 3 were used to compute the 
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static pressure at station 3 by means of figure 2. This pressure 
and the temperature Tj^ are used to compute the density p ex by 

means of the general gas law. 

Computation of C p. r . . - By the method given in aopendix A , 

\~\ 

. ma y be computed from the cooling-surface temperature dis- 

%, - T i t 

tribution and the over -all rate of heat transfer. On the assumption 
that p 2 ^ = P]^* P2^ an( ^ (pV) yP2^. P2^ were calculated and, from 
figure 2, p^, p^, and q 2 were found. From and T^ = , 

the values of p-^, q_,, and p^ were calculated with the aid of fig- 
ure 2. A value for C^. was then computed from equation ($) • 

Computation of exit coefficients. - The exit coefficients C, 
and a could not be separately determined and the value of a^ - C^ 
was consequently calculated from q^, p^, and ^obtained in the 

calculation of and equation til)* From and T^ , figure 2 

was used to find the 7 needed values for and p^. 

Computation of C ^f- by estimation of - C x . - The value of 
^Df . was a l so calculated by the method of reference that is, by 
estimating the exit-loss coefficient and assuming p^ = p^ instead 
of using the measured value of p 7 ^. Equation (17) was then used to 
calculate p^. Also from previous computations p 2 , p 2 , ^2' T! ' and 
M 2 2 were known. The chart (fig. 3) was recalculated for constant 
values of (p 2 - P^)/^ instead of for constant values of Cj)^.. 

new chart was used to find the value of p ! .from (p 2 ~ p-z)/q.2> M 2^' 
T ! ; equation ($) was then used, to compute . 

Computation of Reynolds number. - Reynolds numbers R were 
obtained from the formula 

r.SS ■ 

where 

d hydraulic diameter of fin passage, feet 

[x absolute viscosity of air based on average of average cooling- 
surface temperature and average cooling-air temperature, slugs 
per second per foot 



The 
2 
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Computation of C ^. * - The pressure-loss coefficients Cp^ 

were calculated from the data by means of equation (20), v/hich was 
used in the equivalent form 

From the data it was found that should be calculated from p-, 

°t t 
and Tj^ and that p should equal p lt for the best correlation. 

In the foregoing equation p 0 is the stagnation density. The factor 

( PV) 2 

— * — - — was then calculated and, from the value of this factor, 
P o t P ot 

Po 

— V - " -o 



J~± - M 0 2 was determined from figure 2. 



.... „ _RES.l&IS . Aira_DI3.CUSSI.0N_ l 

Evaluation of Coefficients 

Correlation of heat-transfer coefficients h and U. - The effect 
of altitude on the heat-transf er coefficients h and U is shown in 
figures 8 and 9. The curve for h, the surface heat-transfer coeffi- 
cient, plotted against weight flow of air at the baffle entrance shows 
no effect due to altitude v inasmuch as the curves for U000 and 
ill, 000 feet bracket the spread of the data. The curve for U, the 
■wall heat-transTer coefficient, shows the same characteristics as the 
curve for h. These figures show that the Mach number had no effect 
on the heat transfer for the range covered. An effect of heat load 
was noted in the data but is not shown in figures 8 or 9* The heat- 
transfer coefficients h and U increase slightly with heat-transfer 
load. 

Correlation of pressure-loss function (Ap P Q , r )/r(pV)^. - The 
pressure-loss function (Ap P av )/ "^(pV) is a dimensionless form of 

the conventional pressure-loss function Ap P av /p s « The density p 
is the average of densities upstream and downstream of the cylinder, 
and p s is the air density under standard sea-level conditions. When 

(Ap p av )/|(pV) 2 is plotted against- weight flow pVg as in figure 10, 

the conclusions drawn from the plot of the pressure-loss function may 
bo applied to" the function Ap P av /ps* r ^ v - dimensionless form has . 
tjhe advantage of showing .jcompreGsibliity effects more clearly. 
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The plot of the pressure-loss function Ap P^j^ip^) is shown 

in figure 10 for runs with and without heat transfer. These curves 
shove disagreement for different altitudes and for a variation in 
heat transfer. The rise shown on the high weight-flow region of 
each curve indicates an additional pressure loss caused by compres- 
sibility effects at high Mach numbers. The high Mach numbers occur 
at lower weight flows as the density, altitude increases. Any cor- 
relation method must make it possible to correct for these pressure 
rises and the separation of the curves for different altitudes and 
rates of heat transfer. 

Correlation of pressure-loss function (Ap p P ^|(pV) 2 * - Some 
investigators have proposed correlating the pressure loss 'by means 
of the function Ap p ex , where p ex is the air density at the exit 
of the cylinder. This method was unsuccessfully tried with the data 
of the present tests. In most cases where this method would be 
applied, the baffle exit would have no tailpiece and consequently no 
pressure recovery at the baffle exit. The exit density would then 
be computed from the static pressure at the baffle exit and the stag- 
nation temperature at the tailpiece. The density thus computed is 
lower than any density existing in any part of the test rig used for 
the present tests. The pressure-loss function (Ap P ex )^(pV) 2 was 

computed by means of the density p ex and is plotted in figure 11 
for tests with and without heat transfer. Because p ex decreases 
as compared with p av for increasing pressure losses, the correlation 

of the pressure-loss function (Ap P« x )/|(pV) 2 is better than that 
of (Ap p av )^'i(pV) 2 . Although the correlation of the runs with heat 

transfer (fig. 11(b)) is good, the runs without heat transfer 
(fig. 11(a)) show compressibility effects in that the data for the 
different density altitudes form separate curves. 

Correlation of drag coefficients and exit-recovery coefficients; - 
The drag coefficients Cu f . computed from the test data by the method 

previously given are plotted in figure 12 against the weight flow per 
unit area. The utter lack of correlation between the results for the 
tests without heat transfer (fig. 12(a)) and those with heat transfer 
(fig. 12(b)) or between results at different altitudes is evident. 
The same result is to be observed in figure 13 in 'which the exit- 
recovery coefficient - has been plotted against the cooling- 
air weight flow. 

In order to find the source of error in these results, the 
assumptions involved in this method were examined. With regard to 
the assumption that the total pressure remains unchanged between 
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stations 1 and 2, the data indicated that the pressure drop is so 
small as to be negligible; at high flows the pressure drop amounts 
to about 1 percent of the total-pressure loss* 

The mass flow was calculated from the total and static pres- 
sures at stations 2 and 3 1 to test the validity of the assumption 
that the flow is uniform across each section (one-dimensionalflow) • 
The calculation was made from the factor (pV) /p t p^, which was 
obtained from the pressure ratio p/p^ and figure 2. The ratios 
of the v/eight flows to the weight flows calculated from measurements 
at the orifice are plotted for stations 2 and 3 1 in figures ±h and 
l£, respectively, against the 'weight flows obtained from the orifice 
measurements. At station 2, figure lU shows that the ratio is near 
enough unity for application of the assumption of uniform flow. At 
station 3 1 , however, figure 15 indicates that the measurements are 
not accurate at low flows. In the calculation of any coefficient 
based on measurements at one point at station 3 1 * correlation will 
probably not be obtained. It is partly for this reason that Gp vf ., 
as plotted in figure 12, shows such lack of correlation. 

In order to eliminate the use -of .these inaccurate measurements, 
%)f . ~" was recalculated by the-use-of— arras same d-^-xirt— lros-s— in-s-tead--of 

the measured, loss. The value of a^ compensates to some extent for 
the fact that 6* £ 90°. .If a^ + 2 12. sin 6? is assumed constant 

for various values of G^, then for 6^ = 90 separation from the 
cylinder rear can be expected and no pressure recovery due to' pres- 
sure gradients along. the wall will result. In that case a 7 = 0. 
Then, 



l 3 

or 



A, A, 
az + 2 - sin 9 T = 2 J 



l 3 - 2 -2 (1 - sin 9 3 ) (2$) 



If it is further assumed that 0^=0 and = p^ and if 
= q^ (Ajj/Aj) is substituted in equation (17), then 

% 

This equation applies with an incompressible fluid flowing through 
the sudden expansi.on of a straight duct. 
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The method of calculating ^Df j_ when this loss is used has 
been previously given. The results are plotted in figure 16. Much 
better correlation results than was obtained in figure 12, especially 
at low flows. The data for each altitude form a smooth curve. The 
curves separate, however, as soon as the Mach number, which increases 
with altitude for a fixed Reynolds number, becomes appreciable, indi- 
cating that the effect of compressibility has been overestimated by 
baking too high a loss at the exit; that is, a^ is too low. 

In order to determine how much these irregularities in Cr)^ 
and - affect the over-all total -pressure drop, the values of 
^Df i and a 3 ~ c 3 determined from the tests were used to calculate 
the loss in total pressure to be expected with an incompressible 
fluid by assuming that the density remains unchanged in all equations 
for pressure loss. If all losses are added and if p = constant = p 9 , 

Pi - ■ 1 * ° D fi * c 5 - a ? - % e 5 \rj <«« 

The over-all pressure-loss coefficient C pi is plotted against 
the weight flow per unit area in figure 17, which shows that the 
extreme scatter of the data points in figures 12 and 13 has been suf- 
ficiently reduced to allow fairing of curves through the data points. 
The effect of compressibility may be seen from the fact that the 
curves show systematic differences for the different altitudes although 
the ends of the curves do not have the sharp curvature that is evident 

in the plots of (Ap P av )/|(pV) 2 . These differences are much smaller 

than the entire compressibility effect, as may be seen by comparison 
with figure 10, which is already somewhat corrected for compressibility 
effects. This result would seem to indicate that, for correction of 
the compressibility effect, the division of the pressure loss between 
the baffle channel and the baffle exit is unimportant because the pres- 
sure pj was no doubt in error. In order to investigate this 
hypothesis, the over-all pressure-loss coefficient C pi was computed 
from the estimated value for a^ - (equation (2$)) and the corre- 
sponding values for C Df> (fig. 16) . The results of these calculations 

are plotted in figure 18, which shows that the value for cannot 
be estimated in an arbitrary fashion. 

If the values of and - C* from equations ($) and (1?) 

are inserted in equation (26) in order to determine how the pressure- 
loss coefficient C p . depends on the data, there is obtained 
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' 2 "i 


r p 2 - p 5 


+ P 2 


L q 2 



13 



If the density ratios are estimated from the energy equations for 
very low Mach numbers (M = 0) , then 



2 



(l + T«) 



2T ' 



1 + 



T r 



The value of T ! (see equation (8)) is always less than 0.10 and 
the data indicate that the term involving p^ is never greater than 
0.01. At low Mach numbers, therefore, large' variations in estimates 

C n . For high Mach num- 
* i 

p* will definitely 
C T . . Figure 15 shows that the measurements 
are inaccurate at low Mach numbers but much better at high 



of p^ from its true value will not affect 
bers, however, equation (27) indicates that 
affect the result for 
of 



P 3 

Mach numbers; therefore, at high Mach numbers when the data for p^ 
are used, -fair -cor-relation -o£--C p - -may be- expected-. —At- low-Mach -num- 
bers, the inaccurate values for p^ do not affect the result; cor- 
relation may consequently also be expected in this range . 



inaccurate values for 
relation for C 



G 



and 



With 

at low Mach numbers, cor- 



a 5 - C3 

may be expected; whereas, in the high range, 



correlation may be expected only when accurate values for C 



a 3 ~ C 3 



are obtained. 



and 



Figure 12(a) shows that no correlation of C]) f ^ occurs at low 

Mach numbers and fair correlation occurs at high Mach numbers. No 
accurate correlation in any range is obtained for tests with heat 
transfer (fig. 12(b)) but, if scatter is taken about a mean curve for 
the. high Mach numbers, the variations are small compared with the 
value of C pj and therefore do not show up so prominently. The same 

general features of the C Df . plots (fig. 12) may be observed in the 

plots (fig. 13) and the same comments apply. From the C D 



a 3 " C 3 



was 



data (fig. 16) it may be presumed that the estimate of a^ 
incorrect and resulted, in the calculation of faulty values for pj. 
Comparison of C pi of figure 17 with that of figure 18 indicates that, 

in the range of low Mach numbers for all altitudes, the C~ curves 
are the same with the measured and the calculated values or 1 p^, which 
confirms the deduction that, for low values of Mach number, 0 is 



independent of 



In fact, the lowest-altitude curve is approximately 
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the same for both methods of computing Cp^ except at the very- 
highest Mach numbers. The lack of correlation of the data of fig- 
ure 18 at high Mach numbers confirms the deduction that the esti- 
mated values of are incorrect and influence the values for Cp^ 
and that the measured values for high Mach numbers used in calcu- 
lating C n (fig. 17) are accurate* 
iJ i 

The fact that the pressure measurements in the rear of the cyl- 
inder were very inaccurate at the low cooling-air weight flows indi- 
cates a lack of uniformity in the flow, which contradicts the hypothesis 
upon which the flow analysis was built.. . A second inaccuracy is seen 
in the hypothesis by which is calculated from Cq^. The method 

of calculations is based on the assumption that q ^ ipV^ = i LsZ) 

2 2 L P 

increases linearly along the channel; i(q^ + q2), therefore, repre- 
sents a good mean value for q. This assumption is approximately 
true for low Mach numbers but not for high Mach numbers where a much 

more rapid rate of increase is to be expected. Because '^(^ + ^3) ^ s 

higher than the mean effective value of q in the channel, Cj)« 

1 i 

must be lower than the coefficient that may be expected with a truly 
incompressible fluid. The nonunif ormity of flow in the rear (sta- 
tion 3) will also modify the momentum equations for the baffle channel 
and the baffle exit. For these reasons values of and a^ - Cj 

obtained by the assumption of uniform flow and by means of equations (£) 
and (17) will probably not check estimates from measurements of the 
velocity and pressure distribution along the channel and baffle walls. 
Further experiments are required to determine conclusively the effec- 
tiveness of this method. 

Correlation of the coefficient C p^ calculated by empirical 
method. - A plot of the pressure-loss coefficient Cp^ against the 

Reynolds number is shown in figure 19 • Although the correlation 
shown in this plot is satisfactory, some features show that improve- 
ment is still desirable. The points for the lowest density altitude 
without heat transfer are, in general, lower than the average of the 
points. The tendency for this plot to fall off at high Reynolds num- 
bers indicates that the correction for compressibility effects is 
probably too large. This result might possibly be due to the fact 
that the correction factor applied is too large for the pressure losses 
occurring before station 3« This effect does not occur in the data 
for other density altitudes, possibly because the data for the lowest 
densities may be inaccurate in the high-flow region as a result of 
the unsteadiness of the driving engines. 
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In order to determine the magnitude of these deviations on the- 
actual pressures, the corrected pressure drops were plotted in fig- 
ure 20 in the customary manner. The discrepancies noted in fig- 
ure 19 show little effect on the pressure-loss plot. The shape of 
the plots of figure 19 indicates that a linear curve of pressure- 
loss data is not to be expected. 



General Discussion 

Comparison of the correlation methods* - For any set of pressure- 
loss data, the separation of the curves for different density alti- 
tudes is an effect of compressibility that has not been corrected. 
The effect of compressibility is accentuated at the highest weight 
flow for each density-altitude curve. The effectiveness of each 
method of pressure-drop correlation is characterized by the magnitude 
of the spread of the curves at the four weight flows that are the 
maximums for the different altitude curves. At the end points of 
the curves for U000-, llt.OQO-, 2U,000-, and 32,000-feet density 
altitudes, there are, respectively, two, four, six, and eight points 
on all the curves for comparison. The average deviation from the , 
mean of all the curves at these points is given in the following 
tables in percentage of the mean value of the parameter at that 
point: 



Weight flow 
(lb/(sec) 

(sq ft)) 


Mean deviation, percent 


Method 


Ap P av /|(pV) 2 


Ap P ex /|(P V ) 2 


Analytical 


1U.U 


13 


7.3 


7.3 


21.5 


13 


3.9 


U.7 


29.0 


11 


2.8 


5.0 


40.0 


13 


U.9 


2.1 



Reynolds 


Mean dev.i.ation, 


number, 


empirical method 


R 


(percent) 


O.lUxlO* 


U.8 


.21 


2.U 


.30 


2.1 


.Uo 


1.5 
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These tables show that the use of the variable (Ap p av )/-(pV) is 
least effective, that the use of the variable (Ap P ex )/|(P v ) 2 and 

the analytical method are better and about equally effective, and 
that the empirical method using the Prandtl-Glauert factor is most 
effective and of satisfactory accuracy. 

The involved procedure needed for use of the analytical method 
is a definite drawback. Furthermore, certain information is needed 
for pressure-loss predictions that is not necessary when the empiri- 
cal method is used; the relative amount of heat picked up by the air 
in front of the baffle entrance must be obtained and pressures at 
the baffle-channel expansion point are required in order to evaluate 
separately the coefficients C Df . and - in the pressure-loss 

predictions. The testing technique for determining reliable values 
of p^ is yet to be developed. Computations made for a modern 
engine at an altitude of 10,000 feet show that very high Mach numbers 
occur and that separate values of . and aj - are needed. 

The method is somewhat simplified, however, when no baffle tailpiece 
is used, which eliminates pressure recovery at the baffle exit. In 
that case p^ = pj^, and. pressure measurements are therefore unneces- 
sary in the baffle exit. 

The same drawback applies to the use of the pressure-loss function 

(Ap P ex )/|(pV) 2 ; that is, a measurement of the pressure in the baffle 

exit is required to determine p ex unless a baffle is used that gives 
complete loss of the kinetic energy at the baffle exit. This method 
of correcting the pressure loss for compressibility effects has a 
rational basis. If the entire loss is assumed to occur at the baffle 
exit, if no pressure recovery exists in the baffle tailpiece, and if 
the air temperature in the baffle exit is assumed equal to that behind 

the cylinder, (Ap Peyj/|(p v ) 2 should very nearly account for all 

compressibility effects by equation (17) • If not all the loss in 

energy occurs at the baffle exit and if some pressure recovery is 

present there, this method should give good results if most of the 

pressure loss occurs at the baffle exit and if the rest varies in 
much the same manner as the baffle-exit loss. 

The empirical method of computation of the pressure loss has the 
advantage of being quite simple. There is doubt as to the types of 
flow apparatus to which the empirical method can be applied because 
no rational basis for it exists at present. Test data from engines, 
engine cylinders, and radiators are required before an estimate of 
the reliability of this method can be made. The same confirmation of 
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the reliability of the analytical method is needed, however, because 
there is some doubt that it correctly represented the actual flow 
conditions in even this comparatively simple type of flow path. * 

In order to illustrate the application of the empirical method, 
an example is- computed and plotted in figure 21. For comparison, 
the pressure loss is also computed by means of the function 

(Ap P av y|(P v )^* The barrel temperature is assumed constant at 3£°° F| 

the rate of heat transfer is assumed constant and such as to require 
an air -weight flow vat sea level of 18. $ pounds per second per square 
foot of free-flow area of the barrel-baffle channel. The atmospheric 
pressures and temperatures are modified for isentropic compression 
corresponding to a flight speed of 200 miles per hour. The subsequent 
steps are: 

It Compute cooling-temperature differential - Ti t 

2. Compute the heat-transfer coefficient U (equation \{2h)) 

3. Determine pVg (fig. 9) i compute the exit temperature 

U. Compute the Reynolds number, and determine C pi (fig* 19) 

$. Compute the abscissa of figure 2 and read the compressibility- 



6. Compute the pressure loss from the data of steps 3, k$ and S> 
When the pressure loss is computed by means of the function 



means of the plot of lowest density altitude (fig. 10(b)). An esti- 
mate of Ap is made for determining p av ; the resulting value of Ap 
is used to refine the values for p av and Ap. For less extreme 
cooling conditions, the difference between the two methods is not so 
great. 

Effect of viscosity . - For the general correlation of pressure- 
loss and heat-transfer data, dimensionloss parameters should be used. 
The variable h, for example, should be plotted in terms of the 

Stanton number — ~- against the Reynolds number. In turbulent flow 

°p pV 

the Stanton number varies approximately as the viscosity to the 
l/5 power. For practical application of heat-transfer data to 
altitude-cooling predictions, therefore, the small variation in the 



correction factor 




given. The weight flow is used to determine 
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viscosity caused by the variations in temperature of the cooling-air 
surface film will affect the heat transfer only slightly. Good cor- 
relation of heat-transfer data of a given engine or cylinder can be 
expected, therefore, if the weight flow is used as the correlation 
variable. The same reasoning can be applied to pressure-loss data. 
Because the present pressure-loss data were obtained with and without 
heat transfer, considerable difference should be expected in the film 
viscosity. Consequently, the Reynolds number was used to correlate 
the pressure-loss data; the viscosity of the surface film gave better 
correlation than the cooling-air viscosity in the computation of the 
Reynolds number. The pressure-loss coefficients were plotted against 
the weight flow because use of the Reynolds number effected no increase 
in correlation or significant change in the plot. 

Velocity distribution. - Some improvement might be made in the 
analysis of flow in radiators, for which the assumption of uniform 
flow at the radiator-tube exit is satisfactory, if allowance were made 
in the momentum equation for the nonuniform! ty of flow to be expected 
with a turbulent velocity distribution in the tube. The computed 
kinetic-energy loss at the exit would thus be increased. 



•SUMMARY OF RESULTS 

Based on an analysis of pressure-drop requirements and on experi- 
ments to determine the effect of air compressibility on cooling and 
pressure loss of a baffled cylinder barrel, the following results 
were obtained; 

1. The pressure loss from the front of the cylinder to the baffle 
entrance was very small, as might be expected for any baffled engine 
cylinder. 

2. The method of analyzing the flow processes around a baffled 
cylinder based on the assumption of uniform flow corrected for most 
of the effect of compressibility. 

3. The assumption of uniform flow used in the analytical method 
of predicting pressure losses was not verified by computations from 
pressure measurements between the fins at low flows. 

k* Prediction of the pressure losses across a cylinder at high 
Mach numbers using the analytical method based on uniform flow 
requires knowledge of the exact values of the friction coefficient 
between the fins and the coefficient at the exit of the cylinder. 
At low Mach numbers the division of the pressure coefficients has 
little effect on predicted pressure loss. 
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The use of a fictitious exit density p ex in correcting 
pressure-loss data gives as accurate correlation as the analytical 
method based on the assumption of uniform flow. 

6. An empirical method has been found which gives satisfactory 
correlation (mean deviation of pressure-function curves of about 
3 percent) of the test data on pressure loss of the present cylinder 
and which permits estimates to be made of compressibility and heating 
effects more simply than either the analytical method or the method 
using the baffle-exit density in baffle-flow systems with a tailpiece. 



CONCLUSION 

If engine-cooling equations are based on the cooling-air weight 
flow instead of on the pressure loss, good correlation of data and 
prediction of performance may be expected because test data and 
theoretical estimates showed no appreciable effect of air compres- 
sibility on heat-transfer coefficients in .the ordinary range of engine 
operating conditions* 



Aircraft Engine Research Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio* 
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APPENDIX A 

ESTIMATION OF AIR-TEMPERATURE DISTRIBUTION 

AROUND A BAFFLED CYLINDER BARREL 

In order to determine the pressure loss over the cylinder, it is 
necessary to know the rate of heat transfer to the cooling air between 
stations 1 and 2 H-^ and the rate of heat transfer to the cooling 
air in the baffled section I^. The following analysis presents one 
method of estimating Hj_ and R2. 

For an element of heat-transmitting surface dS, the amount of 
heat transferred per unit time is 

dH = h x (T f - T t ) d5 (Al) 

where 

H rate of heat transfer from front of cylinder to point being 
considered, Btu/sec 

h x local surface heat-transfer coefficient, Btu/(sec) (°F) (sq ft) 

Tf cooling-surface temperature average at any local flow sec- 
tion, °F absolute 

T t local stagnation-air temperature, °F absolute 

The viscous shearing stress slows down the air in the immediate 
vicinity of the wall and increases the temperature there until T, 
the temperature upon which the heat transfer depends, is more prop- 
erly the stagnation temperature than the static temperature (refer- 
ence 11). If W L is the local weight flow between the fins (lb/sec), 
the energy equation is 

H a Cp W L (T t ~ T lt ) (A2) 

In the unbaffled section of the cylinder, W L increases from 
a value of zero at the front stagnation point to the value W at 
the baffle entry. The term T t is eliminated from equations (A2) 
and (Al) and the quantities 
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and 

t = Tf - Ti t = cooling-temperature differential, °F 
are used to obtain 

dH = C (c p W L t - H) d(S/S T ) 

where 

S*£ entire surface for heat transfer, sq ft 

3 surface for heat transfer from front of cylinder to point con- 
sidered, sq ft 

The solution of this differential equation is 



H = exp 



r S/Srp 

Cd(S/S T ) 



Constant + 



,3/S T 

J Cc p Wr t exp 



J Cd(S/S T ) 



d(S/S T )> 
(A3) > 



Because S/Sj = 0 and H = 0 at the front of the cylinder, the 
constant of integration is zero. At the baffle entry, S/Sm is 
equal to a, the ratio of the unbaffled cooling surface of the 
cylinder to the total cooling surface, and H = K^, the rate of 

heat exchange in the unbaffled portion of the cylinder. At the 
baffle exit, S/St = 1> and H = + #2* the rate of heat exchange 
over the entire cylinder. When these boundary conditions are 
inserted in equation (A3), 



Hi ~ exp 



r a 1 r a 

- J Cd(5/S T ) I Cc p W L t exp 



S/S T 



Cd(S/S T ) 



d(S/S T ) (AU) 



Cd(S/ST)j J Cc Y/ L t exp J Cd(S/S T ) 



-S/St 



d(S/S T ) (A*) 



The rest of this appendix derives expressions for H]_ and 
**l' + **2> which are more convenient than equations (Ah) and (A 5) for 
computation from the test data. The quantity C varies little over 
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' / o S/3T Cd(3/5 T ) 
the cylinder except at the very front. The factor e 

! 3Ar Gd(s/S T ) 

and the integral of CW L te'° d(s/s T ) are only slightly 

affected by this variation. Consequently, the quantity C is con- 
sidered a constant in all the integrations. The values found for 
H x and H x + H 2 from the energy equations 

Hi - Wc p :T 2t - T lt ) 

and 

Hi + H 2 = Wc p (T^ - T lt ) 
are substituted in equations (AU) and (A5) to obtain 



T 2 . - Tt. - C e 



-Co f a ^ te C(3/S T ) d(s/3 (A6) 



-c P \ te C(s/s T ) + p J te c(s/s T ) d(s/sT) 



o v a 



(A7) 



T t ,-T tl -e- C <l-> - T lt ) ♦ Ce- C £ 5 te^' d(( ^, (A8, 

If an estimate is made of the weight-flow variation W L /W and 
measurements of t are obtained, G can be determined from the 
over-all temperature rise Tj^ - T-^ and equation (A?) and its 

value used in equation (A6) to find - 

In order to estimate the vreight-flow variation in the unbaffled 
portion of the cylinder, the flow in the front is computed by means 
of the equation for potential flow over a circular cylinder. The 
velocity over such a cylinder is 

V - 2 V 0 sin n (S/S T ) 

where V 0 is the velocity of the undisturbed flow, which reaches a 
maximum at s/s«p = l/2. In the case of the baffled cylinder the 
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maximum value is at S/S? « a. The velocity function is therefore 
modified to take this factor into account 

V = 2 V c sin (n S/2g S T ) 

The constriction at the baffle entrance (S/S-j = a) of the 
flow area to l/n its value for an unbaffled cylinder increases 
the maximum velocity by the factor n, but the velocity in front 
of the cylinder is unchanged. The velocity distribution must 
therefore be multiplied by an even function of S/St* which 
increases from a value of 1 at 5/o«p =0 to a value of n at 

S/S T « c. The assumed function is n |l ~ n - 1 cos (n S/2a St)J 
which gives an approximate weight-flow variation of 

~ = sin (it S/2a S T ) - 2^ sin 2 (s/S T ) (A9) 

From the baffle dimensions, n * For the unbaffled portion of 
the cylinder, a linear cooling-surface temperature distribution is 
assumed ~~ — ~~ " ~~ 

t/U-^-QLZ^P}^ (AID) 

where is the value of t at S/S^ = 0, and t2 is the 

value of t at S/Sy = a. Equations (A9) and (A10) are substituted 

in equation (A6) to obtain 

■ ( T 2 t - - f i + (Srv f2 (c) (A11) 



where 




/ 



NACA ARR No. Eb.G20 



31 



f 2 (C) 5 



a 



'\2a; J 



ft 



ic 2 



Tf 

••2a/ 



n - 1 
2n 



• 2-i 2 
a ! C + la/ i 



2C - + e- Ca 
a 



An. approximation must now be found for the integral involved 
in equation (A8). In the baffled section Wj^'W = 1 and t varies 

very little 5 an extreme case shows 9-percent average deviation from 
the mean. The unbaffled cylinder rear shows a much larger devia- 
tion but the decrease in the effective cooling-air velocity reduces 
the variation of the function (W L /W) t. The value of W L /W is 
therefore taken to be unity and an average value used for t so that 



w L c(s/st) , , x 1 f r n a '\ 



(A12) 



• a 



where 



: 7 is the average cooling-surface temperature of the cylin- 



der from the baffle entrance to the rear minus the inlet -air tem- 
perature. This result is substituted in equation (A8) to obtain 



if 



i t - T it) - h ~ *S - *0 



e C(l-a)_ x 



(A13) 



equation (All) can be plotted in the form =T 2 ^ - T lt.</ *a is a 

function of : (e 0 ^ 1 ""^ - 1 ) and (t x - t 2 )/t a j. This plot is 
shown in figure 22. Equation (A13), hoy/ever, gives a solution for 



the variable 'T?^ - T lf j/ t a 



^ ( 1— <x } 

as a linear function of e° l " ' - 1 



with a y intercept of (T^ - T l4 .j/t a and a slope of 



32 



NACA ARR Ho. EUG20 



Jtj - h]^ - T^] Jt a . From data on the temperature of the cylin- 
der y the straight line (equation (Alj)) can be located on figure 22 
and the intersection point can be obtained y&th the curve of equa- 
tion (All) for the proper value of (ti - t 2 )/t a ; thus, 

( T 2 t ~ T lty* t a can be determined. 
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appendix b 
symbols 



a -f- " P?) ^ 



A cross-sectional area at a station indicated by a subscript 
(except with subscript b) , sq ft 

A b area of cylinder wall at base of fins, sq ft 

Cp specific heat of air at constant pressure, Btu/(lb) (°F) 



C = 



(h x S<p) 



(c p W L ) 

coefficient for friction drag in baffle exit 

coefficient for skin-friction drag from station 2 to station 3 
based on dynamic pressure at station 2 

Cj) f . coefficient for drag from station 2 to station 3 based on 
1 dynamic pressure averaged in baffle channel 

Cp^ pressure-loss coefficient with an incompressible fluid 

C- pressure-loss coefficient with a compressible fluid 

g acceleration of gravity, ft/sec 2 

h^ local surface heat-transfer coefficient, Btu/(sec) (°F) (sq ft) 

h average surface heat-transfer coefficient, based on inlet-air 
temperature, Btu/(sec) (°F) (sq ft) 

H rate of heat transfer from cylinder area between front and 
any point considered, Btu/sec 

rate °f heat transfer to air from front of cylinder to baffle 
entry, Btu/sec 

H 2 rate of heat transfer from baffle entry to baffle exit 
(station 2 to station U), Btu/sec 

M Mach number at a station indicated by subscript (except M a , Mb) 



3U 
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\^ (1 + p') + (p« - 1)] + ^1 (1 - p- 2) 



p static or stagnation (indicated by subscript t) pressure at 
any point in fluid indicated by a subscript, lb/sq ft 

Pj static pressure at any point in fluid with incompressible flow, 
lb/sq ft 

q dynamic pressure at any point in fluid indicated by subscript, 
lb/sq ft 

R Reynolds number 

S surf ace _for_. heat transfer _from..f.r.ont. of : cylinder to point -con- 
sidered, sq ft 

S T total heat-transfer surface of cylinder, sq ft 

S x projected area of cylinder-wail surface or curved part of baffle 
surface back of station 3 on plane of A^, sq ft 

t 5 T f -T lt 



taE \ (t 2 + t x ) 

average value of t between stations 2 and U, °F 

T true air-stream or stagnation (indicated by subscript t) tem- 
perature at any point in fluid indicated by subscript (except 
T b , Tf, and T f , °F absolute 

average temperature of cylinder wall at base of fins, 
°F absolute 

cooling-surface temperature averaged at any local flow section, 
°F absolute 

Tf av average temperature of cooling surface of cylinder, °F absolute 
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T ■ = H2/c p W 2 

U wall heat-transfer coefficient of cylinder, Btu/(sec)(sq ft) (°F) 
V velocity of air at any station indicated by subscript, ft/sec 
W weight of air flowing through baffle, lb/sec 



a ratio of unbaffled cooling surface of cylinder to total cooling 
surface 

Y ratio of specific heats for air (1.39U7) 

Q-z angle between radii of cylinder to cylinder rear and to station 
deg (See fig. 1.) 

[i absolute viscosity, slugs/(sec) (ft) 

p local true or stagnation (indicated by subscript t) density at 
any point indicated by subscript (except p s , p av , and P ex ) * 
slugs/cu ft 

p s standard density of air at 29.92 in. Hg and 60° F, slugs/cu ft 
p av average of densities at stations 1 and h, slugs/cu ft 
p ex fictitious exit density 

P' E P2/P3 

Ap total-pressure drop frcm front to rear of cylinder, 



Ap^ loss in total pressure from front to rear of cylinder with an 
incompressible fluid, lb/sq ft or in. water 

Apj_ loss in total pressure from front to rear of cylinder with an 
incompressible fluid under standard density conditions, 
lb/sq ft or in. water 

Subscripts applicable to A, p, T, V, p, M, and q: 

1* 2, 3> 3 ! , u stations indicated in figure 1 

t. . stagnation condition of gas (absence of subscript t indicates 
true-stream condition) 

0 condition characteristic of entire flow 



f L 



local air weight flow between fins, lb/sec 




36 



NACA ARR No. EUG20 



REFERENCES 

1. Pinkel, Benjamin: Heat-Transfer Processes in Air-Cooled Engine 

Cylinders. NACA Rep. No. 612, 1938. 

2. Pinkel, Benjamin, and Ellerbrock, Herman H. , Jr.: Correlation 

of Cooling Data from an Air-Cooled Cylinder and Several Multi- 
cylinder Engines. NACA Rep. No. 683, 19U0. 

3. Schey, Oscar W., and Pinkel, Benjamin: Cooling of a Radial Engine 

in Flight. Jour. Aero. Sci., vol. U, no. 11, Sept. 1937, 
pp. UU3-U52. 

U. Corson, Blake W., Jr., and McLellan, Charles H. : Cooling Charac- 
teristics of a Pratt & Whitney R-2800 Engine Installed in an 
NACA Short-Nose High-Inlet-Velocity Cowling. NACA ACR No. LljF06, 
19UU. 

$. Becker, John V., and Baals, Donald D. : The Aerodynamic Effects of 
Heat and Compressibility in the Internal Flow Systems of Aircraft. 
NACA _ACR, Sept. 1042. 

6. Meredith, F. ¥.: Note on the Cooling of Aircraft Engines with 

Special Reference to Ethylene Glycol Radiators Enclosed in Ducts. 
R. & M. No. 1683, British A.R.C., 1936. 

?. Winter, H. : Contribution to the Theory of the Heated Duct Radi- 
ator. NACA TM No. 893, 1939. 

8. Weise, A.: The Conversion of Energy in a Radiator. NACA TM 

No. 869, 1938. 

9. Brevoort, Maurice J., Joyner, Upshur T. , and Wood, George P.: 

The Effect of Altitude on Cooling. NACA. ARR, March li'43. 

10. von Karaian, Th. : Compressibility Effects in Aerodynamics. Jour. 

Aero. Sci., vol. 8, no. 9, July 19U-* pp. 337-5$6. 

11. Frankl, F.: Heat Transfer in the Turbulent Boundary Layer of a 

Compressible Gas at High Speeds; Frankl, F., and Voishel, V.: 
Friction in the Turbulent Boundary Layer of a Compressible Gas 
at High Speeds. MAC A TM No. 1032, 19h2. 

12. Emmons, H. W. , and Brainerd, J. G. : Temperature Effects in a 

Laminar Compressible-Fluid Boundary Layer Along a Flat Plate. 
Jour. Appl.'lech., vol. 8, no. 3, Sept. 19U1, PP- A-IO^-A-.UO. 



NACA ARR No. E4620 

s 



Table I 



I 



HI 



jfr giljpfJiiliIiiilEiJ.gll|illl§iiiillil§ ISllgjliiilSSl m». 



S8S 











1 



liiiiisiiiiiii^iiii onmawHw siiiiiEiiiiiii mm 



mmmnmmi ymmmim mmmm mm 



i 



i! 

S3 

!! 
11 



I 

i 



aaaaaaaaaa aaaaaa eaa aaaaaaa mm 

2-1 ! ! i-L 



jjpaaaaaaa aaaaaaaaw aaaaaaa wan 



mm^mW^ 333S353533S3S3333 5S52 

wwwSwwNww«SSSwSS8wgS 222222222232222252 



3333553 



i 

£2 

i 



I 

i 



I 
i 



I 
ft 

a 

r. 



is 

'A 



*> 



I 

5 



8§§888&sssssss8S88SS sssssaassssssassss iiiiiiiiiiiiii 



mmmmm wm®m msm mm 



. 4 S2aSM3aSK8SS8C8 
888888S&8SSS8888S8S3 



2*S222S£i^22*22222S 2SJ3"2^!**l32i322^4 •-<«Hr-i C)(DOO> 



sliiSHBSffil Ills 



II lllliiia 

IV W VV CM W OJ VV SW W W 



2 32233333353353 IsslSSl 

g s2asaaa3333aas sssssss 



SB88SaS8S8fsa.S.^8.3.^. 883S833.^3.:.Z8SS ; ^ SSSj^SSrs.i&ijij 



mmwm mmmmm msaa mm 



iliiiiiiiiisiiiiiiii llllSillllililil lllIIIIIBil iilliil 



iiiiiiiiiiiiiiiiiiSi 



? 



§ 

a" 



NACA ARR No. EYG20 



Table 



TABLE 2. - SUMMARY OP DATA FOR TESTS WITH HEAT TRANSFER 



COMMITTEE FOR AERONAUTICS f*he symbols used are defined In appendix Bj 



Density 
altitude 

tw 


(In. Hg 
abe. ) 


T lt 
°F abs.) 


VPg 
( lb/( see ) 
(sq ft)) 


Pit ' p 4 t 
(In. Hg) 


T 4 t 
(°F abs.) 


P 2t 
(In. Hg 
abs.) 


P 2 
(In. Hg 
abs.) 


P 3't 
(In. Hg 
abs. ) 


P3 

(In. Hg 
abs.) 


P4 
(In. Hg 
abs.) 


?! 


"av **t 
(°F) 


T 


b " T l t 
(°F) 


4,000 


27.80 
28.00 
28.04 
28.02 
27.99 
27.94 
28.15 
28.09 
28.01 


550.6 
551.1 
552 .1 
553.1 
553.6 
555.6 
550.1 
547.1 
553.8 


12.22 
14.97 
18.16 
21.52 
26 .30 
29.55 
36.12 
39.62 
40.89 




i 0.64 
! .70 
1 .01 
1*39 
2.12 
2.75 
4.69 
6.76 
7.60 


630.6 
628.6 
623.1 
620.6 
613.6 
611.6 
597.1 
593.6 
594 .6 


27.77 
27.97 
28.03 
28*01 
27.95 
27.93 
28.13 
28.02' 
27.96 


27.27 
27.24 

26 ".43 
25.56 
24.85 
23.50 
22.13 
21.69 


27.59 
27.75 

27*64 
27.42 
27.21 
27.06 
26.68 
26.79 




27.05 
26.94 
26.50 
25.85 
24 .65 
23.46 
20.85 
17.57 
15.86 


27.02 
26.96 
26.49 
25.90 
24.77 
23.62 
21.13 
16.08 
16.45 




310.56 
302.17 
283.50 
263*39 
250.09 
238.69 
202.06 
200.30 
180.16 


367.00 
365.55 
351.15 
340.09 
321 .75 
310.63 
289.25 
273.38 
246.21 


14,000 


20.47 
20.47 
20.47 
20.48 

20*46 
20.61 
20.53 


540.4 
541.0 
541.3 
542.1 

542 ll 
542.1 
549.2 


8.55 
11.37 
14.54 
17.26 
19.68 
22.31 
26.87 
29.05 




.39 

.61 
.85 
1 1.22 
; 1.64 
! 2.18 
3.44 
4.61 


615.8 
606.2 
597.4 
593.2 
589.1 
585.6 
580.9 


20.48 
20.44 
20*44 
20.44 
20.37 
20.45 
20.60 


20.07 

19.78 ' 

19.39 

18.94 

18.38 

17.82 

16.63 

16.16 


20.38 
20.30 
20.24 
20.12 
19.90 
19.34 
19.74 
19.47 




20.61 
20.28 
19.64 
19.32 
18.63 
17.42 
16.12 
14.18 




19.90 
19.52 
19.06 
16.51 
17.79 
16.91- 
15.14 
13.05 




234.96 
217.77 
198.07 
105. 02 
172. 66 
166.00 
151.53 
142.12 




270.46 
255.09 
239.09 
227.09 
217.54 
211.17 
197.29 
165.68 




















24,000 
* 


14.98 
14.97 
15.03 
15.06 
14.87 
14.84 
14.84 
14.92 


549.6 
551.1 
551.6 
556.1 
557.6 
557.6 
556.6 
559.6 


8.20 
11.49 
16.41 

21.52 
10.09 
13.96 
15.85 
18.26 




.52 
.89 
1.43 
4.10 
.77 
1.35 
1.80 
2.41 




625.6 
616.1 
606.1 
602.1 
650.6 
632.6 
624.6 
622.1 


15.00 
14.98 
14.99 


14.52 
14.09 
13.44 
12.26 
14.13 
13.49 
13.06 
12.66 


14.80 
14.67 
14.67 
14.35 
14.51 
14 .47 
14.33 
14.30 




14.71 
14.06 
13.25 
9.35 
14.15 
13.34 
12.62 
11.69 




14.24 
13.54 
12.76 
8.75 
13.65 
12.84 
12.18 
11.19 




251 .15 
227.34 
204 . 24 
183.06 
308.47 
263.83 
246.77 
230.82 




289.00 
268.54 
246.92 
227.25 
360.65 
317.38 
300.13 
207 . 00 








14.85 
14 ,e2 

14.81 
14.89 










32,000 


11.18 
11.18 
11.19 
11.18 
11.18 
11.24 


537.6 
539.3 
539.6 
540.4 
542.1 
542.8 


8.55 
9.76 
11.37 
12.22 
14.01 
14.37 




.60 
.92 
1.28 
1.66 
2.22 
2.56 




603.3 
604.3 
600.8 
598.9 
598.2 
598.2 


11.19 
11.16 
11.16 
11.15 
13 .17 
11.25 


10.36 
10.11 
9.75 
9.49 
8.99 
8.93 


11.01 
10.95 
10.92 
10.86 
10.81 
10.90 




12.10 
11.77 
11.36 
10.97 
10.26 
9.96 




10.08 
9.67 
9.17 
8.60 
7.67 
7.23 




208.02 
210.69 
201.45 
197.11 
191.92 
1R9.42 




240.92 
245.92 
237.46 
235.04 
231 .79 
720.96 


Density 
altitude 
(ft) 


U 

(Btu/(hr) 
(sq In.) 
(*F)) 


h 

(Btu/(hr) 
(sq In.) 
<°F)) 


<°F) 




t 




APP av 




APPex 






c "*i 


Cpi 

analyt 
leal 




C P1 
(empir- 
ical 




Reynold 

nun*. 




APl a 


(°F) 


1/2 (PV) 




1/2 (PV) 2 


a 3 ' W 3 


method) 


method ) 


berxlO" 5 


(In. H?0) 


4,000 


0.636 
.758 
.878 
1.020 
1.172 
1.271 
1.402 
1.610 
1.505 


0.0861 
.1052 
• 1246 
.1482 
.1728 
1 .1897. 
.2301 
.2520 
.2359 


365.0 
359.5 
540.0 
323.5 
302.5 
286.0 
258.5 
254.5 
225.5 


0.8055 
.8372 
.8529 
.8624 
.9167 
.9371 
1.0948 
1.0412 
1.0687 


236.6 
224.3 
203.5 
189.8 
169.7 
161.4 
127.8 
126.3 
112.2 


1.019 
.880 
.863 
.846 
.850 
.846 
.948 
1.077 
1.098 


0.931 
.802 
.780 
.749 
.735 
.723 
.747 
.783 
.733 




0.4572 
.4595 
.4118 
.4190 
.4522 
.4632 

. .4386 
.3918 
.4399 


0.2347 
.1612 
.1034 
.1070 
.1459 
.1530 
.0904 
.0724 
.0633 


0.6797 
.6039 
.6738 
.5902 
.5959 
.5920 
.5540 
.5828 
.5256 


0.9299 
.^979 
.7741 
.7361 
.7173 
.6980 
• 6913 
.6852 
.6340 


0.1126 
.1365 
.1696 
.2024 
.2501 
.2618 
.3544 
.3905 
.4047 


5.432 
6.994 
9.987 
13.331 
19.405 
23.808 
35.280 
42.074 
41.339 


14,000 


.571 
.726 
.814 
.927 
1.016 
1.098 
1.261 
1.221 


.0753 
.0975 
.1127 
.1305 
.1467 
.1602 
.1885 
.1851 


270.0 
254.0 
238.0 
222.8 
209.8 
220.2 
190.4 
180.8 


.6444 

.7087 
.7563 
.8416 
.8270 
.9977 
.9863 
1.0039 


183.4 
165.2 
144.0 
129.2 
119.7 
110.8 
99.2 
90.1 


1.130 
.986 
.841 
.857 
.866 
.876 
.931 

1.018 


1.041 
.918 
.773 
.766 
.764 
.732 
.739 
.730 




.1157 
.2830 
.4281 
.4544 
.4740 
.6453 
.4615 
.4015 


.0927 
.1844 
.1574 
.1509 
.1933 
.3176 
.1215 
.1039 


.8792 
.8036 
.6315 
.5987 
.6215 
.5747 
.5622 
• 5246 


1.0368 
.9151 
.7647 
.7565 
.7538 
.7343 
.6877 
.6647 


.0827 
.1111 
.1442 
.1722 
.1983 
.2259 
.2747 
.2969 


2.966 
4.626 
6.322 
8.805 
11.417 
14.296 
19.422 
21.937 


24,000 


.515 
.664 
.812 
1.040 
.622 
.788 
.858 
.950 


.0680 
•0900 
.1126 
.1480 
•0834 
.1086 
.1176 
•1354 


284.0 
262.3 
259.6 
221.1 
353.0 
308.5 
291.0 
276.0 




.5926 
.6645 
.7474 
.7879 
.6912 
.7942 
.8454 
.8985 


201.4 
176.1 
151.1 
129.5 
241.5 
197.2 
181.5 
165.9 


1.160 
1.011 

.903 
1.168 
1.100 

.977 
1.010 

.986 


1.061 
.898 
.761 
.785 
.956 
.844 
.837 
.779 




• .2844 
.4033 
.4203 
.3935 
.4806 
.3866 
.4004 
.4217 


.2994 
.2352 
.1532 
.0710 
.3062 
.1716 
.1629 
.1176 


.9172 
.7341 
.6351 
.5797 
.7998 
.6872 
.6647 
.5981 


1.0596 
.8906 
.7372 
.6980 
.9621 
.0306 
.8156 
.7322 


.0779 
.1107 
.1505 
.2121 
.0923. 
.1308 
.1501 
.1739 


2.787 
4.600 
6.847 
12.641 
3.630 
6.332 
8.015 
9.547 


32,000 


•658 
.616 
.699 
.725 
.810 
.850 


.0741 
.0894 
•0944 
.0992 
.1122 
.1151 


238.0 
242.9 
234.1 
232.6 
227.1 
226.0 




.6290 
.6534 
.6929 
.7097 
.7583 
.7788 


165.5 
162.6 
151.4 
146.0 
140.1 
138.4 


.932 
1.098 
1.060 
1.213 
1.162 
1.240 


.856 
.957 
.945 

1.026 
.961 

1.049 




.3764 
,2402 
.2565 
.1614 
.2224 
.1354 


.1802 
.1825 
.1551 
.1250 
.0747 
.0859 


.7040 
.8445 
.8008 
' .8658 
.7545 
.8527 


.8302 
.9445 
.9161 
.9910 
.8890 
.9600 


.0842 
.0960 
.1126 
.1211 
.1392 
.1429 


2.374 
3.520 
4.632 
5.789 
6.825 
7.752 



NACA ARR- 



No. E1G20 



Fig. I 




/ 



NATIONAL ADVISORY 
COMMITTEE FOR AERONAUTICS 

Figure l.— Finned cylinder, baffles, failpiece, and measuring- station 

locations used in tests. 
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Baffle inlet 




Thermocouple 



Baffle exit 



Station 2 
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p Total-pressure tube 
° Static-pressure tube 



Figure 6.* Location of thermocouples and pressure tubes on test 
cylinder. 
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Figure 7.~ Sketch showing method of installing pressure tubes 
between fins. 
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Figure 8.- Variation of surface heat-transfer coefficient 
with cooling-air weight flow. 
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Figure 9.- Variation of over-all heat-transfer coefficient 
with cooling-oir weight flow. 
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(b) With heat transfer. 

Figure 10. — Voriotlon of pressure- loss function with cooling-air 
weight flow. 
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Figure 1 1 .— Variation of pressure-loss function with cooling-air 
weight flow. 
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Figure 12.— Variation of baffle-channel drag coefficient with cooling- 
air weight flow. Drag coefficient calculated from test data. 
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Figure 14.— Ratio of mass flow measured at station 2 to mass 
flow at station 2 obtained from orifice measurements. Tests 
made without heat transfer. 
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Figure 15.— Ratio of mass flow measured at station 3' to mass flow 
at station z' obtained from orifice measurements. Tests made 
without heat transfer. 
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Figure 19.- Effect of Reynolds number on pressure-loss coefficient 
computed by empirical method. 
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Figure 20. ~ Variation of corrected pressure loss with cooling-air 
Weight flow. 
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